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The present status of pentaquark spectroscopy in lattice QCD is reviewed. This talk also includes a brief 
introduction of pentaquark baryons. 



1. Introduction 

The quantum chromo-dynamics (QCD) may 
not preclude the presence of the multi-quark 
hadrons such as tetraquark (qqqq), pentaquark 
{qqqqq), dibaryon (qqqqqq) and so on, because of 
the color confinement. Especially, we are let to 
be interested in exotic multi-quark hadrons, which 
should have exotic quantum numbers. Now, let 
us address ourself to the pentaquark state. Con- 
sider the SU(3) flavor case. Pentaquark states 
should form six different multiplets: 
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where subscripts in the right hand side denote 
the number of degeneracy in each multiplet. The 
first three multiplets are common in the case of 
usual baryons. However, the last three multiplets; 
antidecuplet, 27-plet and 35-plet, are distinct ir- 
reducible representations since those multiplets 
have an apparent exotic quantum-number such as 
strangeness +1. Needless to say, S — +1 baryon 
can not be accommodated by usual baryons. 
Other possible exotic quantum-numbers are rep- 
resented as stars in Fig. ^ Of course, one cannot 
predict which multiplet is preferred for the possi- 
ble S — +1 pentaquark baryon, within the group 
theoretical argument. The Skyrme model [X] and 
the chiral soliton model [2], however, predict that 
the lowest S = +1 state appears uniquely in the 
antidecuplet and its spin and parity is spin-half 
and positive parity. 

Recently, the LEPS collaboration at Spring-8 
has observed a very sharp peak resonance in the 
K~ missing-mass spectrum of the 771 — > nK + K~ 



reaction on 12 C 3 j. The observed resonance 
should have strangeness +1. Thus, O + (1540) 
cannot be a three quark state and should be 
an exotic baryon state with the minimal quark 
content uudds. The peak position is located at 
1540 MeV with a very narrow width. Those are 
quite consistent with the chiral-soliton model's 
prediction jzj . This discovery is subsequently con- 
firmed by other experiments |4I5| . Experimen- 
tally, spin, parity and isospin are not determined 
yet. Non-existence of a narrow resonance in pK + 
channel indicates that possibility of J = 1 has 
been already ruled out |3I5| 

Many theoretical studies of pentaquarks are 
also triggered by the discovery of the O + (1540). 
I introduce the most reputed model proposed by 
Jaffe and Wilczek In the naive quark mod- 
els, the low-lying pentaquark state should have 
spin-1/2 and negative parity. However, in this 
case, the pentaquark baryon just falls apart into 
KN in a S-wavc. It is difficult to explain its 
very narrow width. Jaffe and Wilczek propose 
a simple idea to flip the parity of the low-lying 
pentaquark. Suppose there is the strong diquark 
correlation. The spin-0, color triplet and flavor 
triplet diquark would be favored within the sim- 
ple one gluon exchange. The pentaquark can be 
composed of two identical bosons (diquarks) and 
one antiquark. However, the anti-symmetrization 
in terms of color, requires relative odd number's 
angular-momentum between the pairs of identi- 
cal bosons. Otherwise, the wave function of the 
pentaquark state should be vanished. Resulting 
parity of the low- lying pentaquark is same as the 

J It should be noted, however, that the experimental ev- 
idence for the 0+(154O) is not very solid yet since there 
are a similar number of negative results to be reported [H]. 




Figure 1. Weight diagrams for possible pentaquark multiplets. Stars represent the states which have 
exotic quantum numbers (electric charge and strangeness). 



chiral soliton model [7]. In addition, the S = +1 
baryon is uniquely assigned to antidecuplet in this 
description. 

Consequently, correlated quark models, i.e. the 
diquark model, may accommodate the positive- 
parity isosinglet pentaquark Q(uudds) as same 
as the chiral soliton model. However, there are 
essential differences between the diquark model 
and the chiral soliton model. In correlated quark 
models, one cannot fully discriminate the other 
multiplet such as octet. If apart from the SU(3) 
flavor limit, the mixing between octet and antide- 
cuplet should be taken into account. Jaffe and 
Wilczek advocate the ideal mixing case, which is 
favored in the vector meson spectrum [7|. The 
ideal mixing provides the different prediction for 
the mass of possible exotic state, which 

should have isospin-3/2. Furthermore, the Roper 
resonance, N(14AQ), can be accommodated in the 
diquark model owing to the ideal mixing 7 . 

A candidate of the isospin-3/2 pentaquark 
state has been reported by the NA49 collabora- 
tion They observed the sharp peak resonance, 
which has exotic quantum numbers, in S-7T miss- 
ing mass analysis. The observed mass is located 
somewhat between predictions made by the chi- 
ral soliton model and the diquark model. In this 
experiment, the isospin partner, S!^ 2 is also ob- 
served. It should be pointed out that this discov- 
ery is not confirmed yet by other experiments [HJ ■ 

I shortly mention about charm or bottom ana- 
log of the 9 state. If anti-strange quark is re- 



placed by a heavier antiquark, what is going to 
be happened. Several models predict that charm 
or bottom analog of the pentaquark is expected to 
be a bound state |7I9| . Recently, the discovery of 
the exotic anti-charmed baryon (uuddc) has been 
announced by the HI collaboration ^H]- How- 
ever, the observed state is higher than the DN 
threshold and also the D*N threshold against 
several model predictions 7 9. Again, this dis- 
covery has not been confirmed yet by other ex- 
periments. 

2. Lattice pentaquark spectroscopy 

If the pentaquark baryons really exist, such 
states must emerge directly from first principles, 
QCD. Of course, what we should do is to confirm 
the presence of the pentaquarks by lattice QCD. 
Experimentally, it is rather difficult to determine 
the parity of the 6+ (1540). Thus, lattice QCD 
has a chance to answer the undetermined quan- 
tum numbers before experimental efforts. Lattice 
QCD has also a feasibility to predict the masses 
for undiscovered pentaquark baryons. I stress 
that there is substantial progress in lattice study 
of excited baryons recently Especially, the 

negative parity nucleon AT* (1535), which lies close 
to the + (154O), has become an established state 
in quenched lattice QCD |11I12| . Here I report 
that quenched lattice QCD is capable of studying 
the 9+(1540) as well. 

Indeed, it is not so easy to deal with the 
qqqqq state rather than usual baryons (qqq) and 
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mesons (qq) in lattice QCD. The qqqqq state 
can be decomposed into a pair of color singlet 
states as qqq and qq, in other words, can de- 
cay into two hadron states even in the quenched 
approximation. For instance, one can start a 
study with a simple minded local operator for the 
+ (154O), which is constructed from the product 
of a neutron operator and a K + operator such 
as O — s abc (d^C'^ 5 u b )d c (s e j 5 u e ). The two-point 
correlation function composed of this operator, 
in general, couples not only to the 8 state (single 
hadron) but also to the two hadron states such as 
an interacting KN system |13I14| . Even worse, 
when the mass of the qqqqq state is higher than 
the threshold of the hadronic two-body system, 
the two-point function should be dominated by 
the two hadron states. Thus, a specific opera- 
tor with as little overlap with the hadronic two- 
body states as possible is desired in order to iden- 
tify the signal of the pentaquark state in lattice 
QCD. Once one can identify the pentaquark sig- 
nal in lattice QCD, to determine the parity of the 
+ (154O) is the most challenging issue at present. 
Thus, it is necessary to project out the parity 
eigenstate from given lattice data precisely. In 
the following subsections, I discuss three related 
issues 

2.1. Estimation of the KN threshold 

The experimentally observed + (154O) state is 
clearly a resonance state. However, its mass is 
near the KN threshold. We could manage to 
calculate the pentaquark as a bound state if its 
parity were positive. 

I recall that all momenta are quantized as 
pL = 2itI (I € Z 3 ) on lattice in finite volume (the 
spacial extent L) with the periodic boundary con- 
dition (PBC) ? Thus, the spectrum of energies of 
two hadron states such as KN states with zero 
total momentum should be discrete and these en- 
ergies are approximately equal to values, which 
are evaluated in the noninteracting case: 

E« N = ^M%+pl + ^Ml+pl (1) 

2 Of course, this (juantization condition may change to 
pL + 2<5 (p) = 2-nl with a scattering phase shift S(p) by 
an interaction between two hadrons. However, the KN 
channel has very weak interactions as is well known exper- 
imentally. One may omit this shift for a crude estimation. 
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Figure 2. The S-wave and P-wave KN threshold 
energies on a lattice of spatial extent L. In the 
case of L < 4.6 fm, the mass of the + (154O) is 
lower than the P-wave KN threshold. 



where p n = i/n • 2tt/L and n 6 Z. The posi- 
tive parity O state decay into KN in a P-wave 
where the KN system should have a nonzero rel- 
ative momentum. The P-wave KN threshold is 
simply estimated at an energy level E\, which 
is evaluated with the smallest nonzero momen- 
tum pi — 2tt/L in Eq. This energy level E\ 
can be lifted by decreasing spatial extent as de- 
picted in Fig. [3 while the lowest energy level Eo, 
which corresponds to the S-wave KN threshold, 
remains unchanged. The level crossing between 
E\ and the mass takes place around 4.6 fm in 
this crude estimation. It implies that the posi- 
tive parity © state may become a bound state in 
the typical lattice size of currently available lattice 
calculations, i.e. L sa 2 — 3 fm. 

2.2. Choice of operators 

For the case of the negative parity state, the 
presence of the KN scattering state must com- 
plicates the study of pentaquarks in lattice QCD. 
One should choose an optimal operator, which 
couples weakly to the KN scattering state, in 
order to access the pentaquark state above the 
(S-wave) KN threshold. 

For this direction, I would like to recall that 



4 



2.0 



1.5 



1.0 



0.5 



0.0 



8 



p=6.0, V=16 3 x32, DWF 
_i i i i 



0.0 0.1 0.2 0.3 0.4 0.5 

(aMJ 2 



Figure 3. Comparison of the fitted mass from 
(Of Of)_(o), (O2O2) (°) and the cross correla- 
tion (o?o$ + 0^6^) (*) ng. 



internal quark fields as 



(2) 



where 77 = ±1. However, due to the relation 
q(+) = 7 5 cK~) for the local baryon operator, the 
resulting two-point correlation functions are also 
related with_each other as (C)W (0)) = 

-7 5 (O(-)(a;)O(-)(0)>75. This means that the 
two-point correlation function composed of the 
local baryon operator can couple to both positive- 
and negative-parity states. However, I note 
that anti-particle contributions of opposite par- 
ity states is propagating forward in time. Thus, 
the +/— parity eigenstate in the forward propa- 
gating contributions is obtained by choosing the 
appropriate parity projection (1 ± 7774)72, which 
is given in reference to the intrinsic parity of op- 
erators, r\. Details of the parity projection are 
described in Ref. |12| . 



the less known observation in the spectroscopy 
of the nucleon. There are two possible inter- 
polating operators for the I = 1/2 and J p = 
1/2+ state; Of = e abc [u T a C l5 d b ]u c and Of = 
£ohc[Uo C , db]75M c , even if one restricts operators 
to contain no derivatives and to belong to the 
(|,0) © (0,i) chiral multiplet under SU(2) L ® 
SU(2)r ^2]- Of course, two operators have the 
same quantum number of the nucleon. The first 
operator Of is utilized conventionally in lattice 
QCD since the second operator Of vanishes in 
the non-relativistic limit. It implies that the sec- 
ond operator is expected to have small overlap 
with the nucleon | (0|Of |Nucl) | w 0. Indeed, the 
mass extracted from the correlator constructed by 
the second operator Of exhibit the different mass 
from the nucleon mass |12I15| as shown in Fig. [3] 
The operator dependence on a overlap with de- 
sired state is evident, at least, in the heavy quark 
regime, while the cross correlation suggests that 
the small overlap with the nucleon might be no 
longer robust in the light quark regime where is 
far from the non-relativistic description |12| . 

2.3. Parity projection 

The intrinsic parity of the local baryon operator 
can be defined by the parity transformation of 



3. First exploratory studies 

3.1. Local pentaquark operators 

As I remarked previously, an optimal opera- 
tor, which couples weakly to the KN scatter- 
ing state, would be required to explore the pen- 
taquark baryons in lattice QCD. For this purpose, 
two types of local pentaquark operator for isos- 
inglet state are proposed in the first two studies. 
One is a color variant of the simple product of 
nucleon and kaon operators, 

o\~J = e abc [u^C^d b ]{u e {s e ^d c )-(u <-> d)}, (3) 

which is proposed by Csikor et al. {Ijg\ 

The other is proposed by Sasaki ^7] as in a 
rather exotic description guided by the diquark- 
diquark-antiquark structure: 

where r^a = 1, 75,757m (but Ti ^ T 2 ) and the 
superscript "77" stands for the intrinsic parity of 
the operator. There are three kinds of diquark- 
diquark-antiquark operator in this description, 
which are useful for the extended study with the 
3x3 correlation matrix analysis. More details 
of construction of the diquark-diquark-antiquark 
operator are described in Ref. [T7| . 
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Figure 4. Masses of the isosinglet S = +1 
baryons with both positive- and negative-parity 
as functions of pion mass squared |16I17) . The 
experimental value for 8 + (1540) is marked with 
a star. 



In an exploratory study one may assume that 
those interpolating operators have smaller over- 
lap with the K A^-scattering state than the simple 
product of nucleon and kaon operators, at least, 
in the heavy quark regime. Because, in the non- 
relativistic limit, all of them give rise to the dif- 
ferent wave function from the KN two-hadron 
system. 

3.2. Results 

The first two lattice studies were performed 
with the Wilson gauge action and the Wilson 
fermion action at the almost same box size L ~ 
2.0 — 2.2 fm. The lattice spacing Csikor et al. use 
is rather coarse than that of Sasaki, but their cal- 
culation was employed with relatively lighter pion 
masses (M„ ~ 0.4—0.6 GeV). The main difference 
between two studies is the choice of pentaquark 
operators. 

After some initial confusion about the parity 
assignment 3 , both calculations agreed that the 
lowest state of the isosinglet S = +1 baryons has 
the negative parity as shown in Fig. ^ The main 
results from the first two exploratory studies can 



3 See a footnote in Ref. 1161 . 



be summarized as follows. 

• The current lattice simulations seem to give 
no indication of a pentaquark in the posi- 
tive parity channel to be identified with the 
6+(1540). 

• The negative parity channel can easily ac- 
commodate a pentaquark with a mass close 
to the experimental value. 

Therefore, both authors conclude that the ex- 
ploratory lattice study favors spin-parity (1/2) - 
and isospin for the O + (1540). In Ref. |TZ], the 
anti-charmed analog of the 6 state was also ex- 
plored. It is found that the Q c (uuddc) lies much 
higher than the DN threshold, in contrast to sev- 
eral model predictions |7I9| . More detailed lat- 
tice study would be desirable to clarify the signif- 
icance of those observations. 

4. Subsequent lattice studies 

There are four subsequent lattice studies of 
pentaquark spectroscopy to be found in the lit- 
erature |18ll9l2()l21j . Other two preliminary re- 
sults had been also reported at some confer- 
ences |22I23| . I give a short review of those results 
as follows. 

Recently, Kentucky group performed their sim- 
ulations near the physical pion mass region with 
overlap fermions ^Hj- However, they choose the 
simple minded operator as the product of nucleon 
and kaon operators to explore the pentaquark 
baryons. Instead, the sequential constrained fit- 
ting method is applied in their analysis to disen- 
tangle the pentaquark signal from towers of KN 
scattering state. They also check carefully the 
volume dependence of the spectral weight. They 
claim as follows: Only towers of KN scatter- 
ing state are seen in the negative parity channel. 
They confirmed that ground state in either par- 
ity channels has a characteristic volume depen- 
dence on the spectral weight, which should have 
1/L 3 dependence for two particles. They also 
confirm that the ghost contribution from KNrj 1 
state appears from pion mass less than around 
0.3 GeV in the positive parity channel. Their 
final conclusion is that there is no sign of pen- 
taquark signal in either parity channels in their 
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Table 1 

Summary of the present status of lattice pentaquark spectroscopy 



author(s) and reference 


signal 


parity of pentaquark 


kind of operator 


Csikor et al. JS| 


Yes 


negative 


color variant of KN 


Sasaki [T7] 


Yes 


negative 


diquark-diquark- antiquark 


Mathur et al. [IB] 


No 


N/A 


simple KN 


Chiu-Hsieh QU 


Yes 


positive 


diquark-diquark- antiquark 


Ishii et al. |2"U] 


No 


N/A 


diquark-diquark- antiquark 


Alexandrou et al. |2] 


Yes 


negative 


diquark-diquark- antiquark 


MIT group E2] 


Yes 


negative 


diquark-diquark- antiquark 


YITP group ESj 


Yes 


negative 


simple KN and color variant of KN 



calculation. However, it seems that their results 
are consistent with the experimental fact that the 
+ (154O) state has not yet been found in the KN 
scattering data as an elastic resonance |24) . 

Another negative results against the first 
two studies are reported by TIT group [2U| . 
First, they trace calculations of Ref. ^7] on an 
anisotropic lattice with the 0(a) improved Wil- 
son fermion action. They confirm that the lowest 
energy state appears in the negative parity chan- 
nel. In their study, a new method is proposed to 
lift up the S-wave KN threshold by imposing a 
hybrid boundary condition (HBC) in the spatial 
direction and it is also applied in their calcula- 
tion. They found that the plateau of the low- 
est energy state is shifted in the effective mass 
plot by changing the spatial boundary condition 
from PBC to HBC as is expected in the case of 
the two hadron states. However, their calcula- 
tions are employed with smeared sources, which 
are apparently optimized for the ground state of 
usual hadrons, i.e. nucleon and kaon. This opti- 
mization may enhance the signal of the KN scat- 
tering states rather than that of the pentaquark 
state. Indeed, the resulting effective mass plots 
for the pentaquark do not show any effectiveness 
of their smearing in the earlier Euclidean time 
region, in contrast to the case of usual baryons. 
However, of course, their method, i.e. the hybrid 
boundary method, is quite promising to distin- 
guish between the pentaquark state and the KN 
scattering state. 

There are several positive results for the first 



two studies. Cyprus group performed careful 
studies concerning the volume dependence of the 
spectral weight and then found that their ob- 
served state seems to be a single hadron state, 
i.e. the pentaquark state |21|. MIT group for- 
mulates eight possible local operators based on 
the diquark-diquark-antiquark structure to per- 
form the 4x4 correlation matrix analysis in both 
isosinglet and isotriplet channels as a extension 
of Sasaki's proposal |22|. Their preliminary re- 
sults are currently obtained from double expo- 
nential fits on data for each single operator. They 
confirmed results of the first two studies at their 
exploratory stage. YITP group has started the 
2x2 correlation matrix analysis with the simple 
minded operator as the product of nucleon and 
kaon operators and its color variant [23]. Their 
preliminary results support that there is an in- 
dication of the presence of the pentaquark state 
near the lowest KN scattering state in the nega- 
tive parity channel. 

Finally, I comment on Chiu-Hsieh's results. In 
the first version of Ref. > they initially misun- 
derstood the parity assignment. Then, their re- 
sults are completely opposite to ones of any other 
studies where the lowest energy state appears in 
the negative parity channel. In the second ver- 
sion, they corrected this apparent error and rean- 
alyzed their data. However, their final conclusion 
remains unchanged. They insist that the spin- 
parity of the 6+(1540) state is most likely (1/2)+. 
This conclusion seems to attribute to their crude 
chiral extrapolation. 
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5. Summary and Outlook 

Table ^ represents a summary of the present 
status for each lattice calculation. The first con- 
clusion of the first two studies as summarized in 
Sec. 13.21 is confirmed by subsequent lattice stud- 
ies. The currently important issue is whether or 
not to establish the presence of the + (154O) in 
the negative parity channel. It is necessary for 
this to disentangle the pentaquark signal from 
the KN scattering states completely . The cor- 
relation matrix analysis is strongly required to 
separate the KN scattering state and isolate the 
pentaquark state. We also should check the vol- 
ume dependence of the spectral weight. Proba- 
bly the hybrid boundary condition is helpful to 
identify the pentaquark state as a single hadron 
state. Needless to say we ought to try non-local 
types of pentaquark operator in order to verify 
whether there is no indication of the O + (1540) in 
the positive parity channel. 

Finally I stress that all present results should 
be regarded as exploratory. Indeed, much detail 
studies are in progress in each group. Thus, I had 
better to conclude that the following questions 
still remain open: 

• Does the spectrum of QCD possess the 
6+(1540)? 

• What is spin and parity of the O + (1540)? 

• Are there other pentaquark baryons, e.g. 
the charm (bottom) pentaquark or the spin- 
orbit partner of the 9+ (1540)? 

There are many exiting issues to be explored. 

Acknowledgments 

This talk has benefited from conversation and 
correspondence with T.-W. Chiu, Z. Fodor, S. 
Katz, F.-X. Lee, K.-F. Liu, N. Mathur and J. 
Negele. The author is supported by JSPS Grant- 
in-Aid for Encouragement of Young Scientists 
(No. 15740137). 

REFERENCES 

1. A. V. Manohar, Nucl. Phys. B 248, 19 (1984); 
M. Chemtob, Nucl. Phys. B 256, 600 (1985). 



2. D. Diakonov, V. Petrov and M. V. Polyakov, 
Z. Phys. A 359, 305 (1997). 

3. LEPS Collaboration, Phys. Rev. Lett. 91, 
012002 (2003). 

4. DIANA Collaboration, Phys. Atom. Nucl. 66, 
1715 (2003); CLAS Collaboration, Phys. Rev. 
Lett. 91, 252001 (2003). 

5. SAPHIR Collaboration, |hep-ex/ 0307083 

6. For a recent experimental review, see 
K. Hicks, hep-ph/0408001, references therein. 

7. R. L. Jaffc and F. Wilczck, Phys. Rev. Lett. 
91, 232003 (2003). 

8. NA49 Collaboration, Phys. Rev. Lett. 92, 
042003 (2004). 

9. F. Stancu, Phys. Rev. D 58, 111501 
(1998); M. Ka rliner and H. J. Lipkin, 
hep-ph/0307343 

10. HI Collaboration, Phys. Lett. B 588. 17 
(2004). 

11. For recent reviews, see S. Sasaki, 
Prog. Theor. Phys. S uppl. 151, 143 
(200 3), JnucT- th/0305014; C. Morningstar, 
| nucl- th/0308026; D. B. Leinweber et ai, 
nucl-th/0406032 and references therein. 

12. S. Sasaki, T. Blum and S. Ohta, Phys. Rev. 
D 65, 074503 (2002). 

13. M. Luscher, Commun. Math. Phys. 105, 153 
(1986). 

14. M. Fukugita et al, Phys. Rev. D 52, 3003 
(1995). 

15. K. C. Bowler et ai, Nucl. Phys. B 240, 213 
(1984); D. B. Leinweber, Phys. Rev. D 51, 
6383 (1995). 

16. F. Csikor, Z. Fodor, S. D. Katz and 
T. G. Kovacs, JHEP 0311, 070 (2003), 
|hep-lat/ 0309090. 

17. S. Sasaki, Phys. Rev. Lett. 93, 152001 (2004), 
|hep-lat/03100T4l 

18. N. Mathur et a/., |hep-ph/0406196| 

19. T. W. Chiu and T. H. Hsieh, hep-ph/0403020. 

20. N. Ishii et al. , |hepd at /0408 0307" 

21. C. Alexandrou et. a/., |hepda"t /0409065 

22. J. Negele, Talk presented at QNP2004. 

23. T.T. Takahashi, Talk presented at Pcn- 
taquark2004. 

24. R. A. Arndt et ai, Phys. Rev. C 68, 042201 
(2003) [Erratum-ibid. C 69, 019901 (2004)]. 



